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We report inelastic neutron scattering measurements of the phonon spectra in a pure powder 
sample of the multiferroic material BiFeOs. A high-temperature range was covered to unravel the 
changes in the phonon dynamics across the Neel (Tjv ~ 650 K) and Curie (Tc ~ 1100 K) temper- 
atures. Experimental results are accompanied by ab-initio lattice dynamical simulations of phonon 
density of states to enable microscopic interpretations of the observed data. The calculations repro- 
duce well the observed vibrational features and provide the partial atomic vibrational components. 
Our results reveal clearly the signature of three different phase transitions both in the diffraction 
patterns and phonon spectra. The phonon modes are found to be most affected by the transition at 
the Tc- The spectroscopic evidence for the existence of a different structural modification just below 
the decomposition limit (Td ~ 1240 K) is unambiguous indicating strong structural changes that 
may be related to oxygen vacancies and concomitant Fe 3+ — > Fc 2+ reduction and spin transition. 

PACS numbers: 75.85.+t, 78.70.Nx, 71.15.Mb, 63.20.dk 



Multiferroism is the combination of two or more of the 
following properties: ferroelectricity, ferromagnetism and 
ferroelasticity. In particular, tuning the electric or mag- 
netic properties by applying a magnetic or an electric 
field, respectively, is at the basis of new applications in 
microelectronics, spintronics, data storage and comput- 
ing hardware^—. This explains the strong interest in mul- 
tiferroic materials from a practical point of view. They 
unfortunately do not abound in nature and are often dif- 
ficult to synthesize technologically. Within the most im- 
portant ABO3 perovskite family the scarcity of multifer- 
roic properties is mainly due to the (i-electrons of the B 
cation, which on one side are necessary for magnetism 
but on the other side tend to reduce the tendency for an 
off-centering ferroelectric distortion^. A remarkable ex- 
ception is bismuth ferrite BiFeOs (BFO). This material 
seems to be unique in exhibiting multiferroic behavior 
above room temperature. It is simultaneously ferroelec- 
tric, antiferromagnetic and ferroelastic, and is subject to 
crystallographic distortions. The mixture of these mul- 
tiple order parameters involved in its phase transitions 
is believed to lead to a very rich phase diagram^. The 
physical behavior of BFO is characterized by two largely 
different transition temperatures; antiferromagnetic with 
T N ~ 600 - 650 K and ferroelectric with T c ~ 1050 - 
1100 K, and a polarization that can reach a very large 3 ^ 



value ~ 100 /zC/m 2 . The large value of the Curie temper- 
ature could be at the origin of the observed colossal spon- 
taneous polarization through a strong ferroelectric order- 
ing!. The polarization in BFO is mainly caused by the Bi 
6s 2 lone pair, so that ferroelectricity is induced primarily 
via this site (strong cation displacements). Another re- 
markable aspect of BFO is the occurrence of FeOe octahe- 
dra tilts accompanying the strong Bi 3+ displacements at 
room temperature^. Below T n the magnetic ordering is 
non-collinear with a G-type-like cycloidal spin-spiral ar- 
rangement and the structure is incommensurate^ with 
an extremely long period of 620 A. The coexistence of 
high-temperature magnetism and strong ferroelectricity 
makes BFO a material with physical properties of prac- 
tical importance^. Due to these exceptional character- 
istics, BFO has become one of the most studied multi- 
ferroic materials leading to a wealth of experimental and 
simulation results covering both structure and dynamics 
— ilir.il . The keen and continuous interest is equally mo- 
tivated by the controversies raised by a priori incompat- 
ible experimental findings: determination and interpre- 
tation of phase diagram, occurence and nature of phase 
transitions, assignment and characterization of frequency 
modes from different spectroscopy techniques (Raman, 
infrared, ...)2 i id i lfi im &. 

Investigations of the vibrational dynamics could help 
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to lift at least some of the controversies. Numerous 
Raman and infrared measurements 3 - ' 13 ' 14 ' 16 ' 20 show that 
the different phase transitions leave clear signatures in 
the vibrational modes. The limitation of optical spec- 
troscopy to the Brillouin zone centre, however, prevents 
obtaining a complete picture of the interplay of dynam- 
ics and structural instabilities. In particular zone bound- 
ary modes could contribute to the unusual properties of 
BFO like the non-centrosymetric positional shift, which 
induces ferroelectricit y 20 ' 21 . The strong debate concern- 
ing the mode assignments from Raman and infrared spec- 
troscopie a 13 ' 16 ' 22 ' 23 just reflects the difficulties in estab- 
lishing the structural phase diagram. It thus seems ob- 
vious that the structural information obtained from dif- 
ferent diffraction techniques should be adequately com- 
plemented by spectroscopic investigations covering the 
full Brillouin zon o 14 ' 24 " — . A concrete argument in fa- 
vor of adopting such techniques, like inelastic neutron 
scattering (INS), arises from the fact that diffraction ex- 
periments allow to capture only the time-averaged struc- 
ture with the consequence of missing subtle structural 
instabilities and distortions that pervoskites are often 
subjected to. At complete odds with the importance of 
BFO, no INS measurements across both Neel and Curie 
transition temperatures have been attempted up to now. 
The present letter aims at closing this gap by probing 
phonon dynamics over the entire Brillouin zone and over 
a wide high-temperature range covering both the Tjv 
and Tc- We devote particular attention to the cross- 
ing of the Curie temperature T^ where the ferroelectric- 
to-paraelectric transition occurs and concomitantly the 
electric polarization should show a sudden drop. Further- 
more, within this high temperature domain some impor- 
tant changes in the electronic properties occur like the ob- 
served metal-to-insulator transition around 1200 K and 
a presumably high-spin to low-spin transition 3 - ' 12 ' 16 ' 24 ' 25 , 
resulting from the reduction of Fe 3+ to Fe 2+ . The lat- 
ter, if it involves a large fraction of the ions, can lead to 
an increase in the cell volume and, therefore, should be 
accompanied by a strong dynamical signature. 

Polycrystalline samples of BiFeOs were prepared by 
solid state reaction of stoichiometric amounts of Bi2 03 
(Alfa, 99.5 %) and Fe 2 3 (Alfa, 99.5 %). Bi 2 3 was 
heated to 800 °C in order to avoid any carbonate or oxy- 
carbonate in the starting material. Fe203 was heated to 
around 300 °C in order to eliminate any adsorbed wa- 
ter molecules. Appropriate amounts of Bi 2 03 and Fe203 
were thoroughly homogenized using agate mortar pes- 
tle with acetone as grinding medium. The well homog- 
enized powder sample was pelletized to a disk of about 
1" diameter and 0.5" thickness and heated to 100-150 
°C to remove any residual organic component. Further 
the pellet was inserted directly into a furnace with a 
preset temperature of 700° C and soaked for about 4 h 
and then quenched to ambient temperature by removing 
the sample from the hot zone of the furnace. At this 
stage the majority of sample formed as rhombohedral 
BiFeOs. Further the pellet was crushed to powder and 



repelletized and heated in a similar manner at 850 °C for 
4 h and quenched to ambient temperature. The sample 
obtained after this heat treatment was characterized by 
XRD recorded on a Philips Xpert-pro diffractometer in 
the anglular range of 10-110° with a step width of 0.02° 
using Ni filtered Cu-K radiation. All the observed reflec- 
tions could be assigned to BiFeOs (rhombohedral R3c). 
The powder XRD data was further analyzed by Rietveld 
refinement with the reported structural data for BiFeOs. 
The final refined unit cell parameters a = 5.5775(2) and 
c = 13.8648(5) A are in good agreement with parameters 
reported earlier for BiFeOs . 

The INS measurements were performed at the Institut 
Laue Langevin (ILL) (Grenoble, France) on 8 grams of 
powder sample prepared as described above. We used 
the cold neutron time-of-flight spectrometer IN6 operat- 
ing with an incident wavelength = 4.14 A (E^=4.77 
meV), which provides a good resolution within the con- 
sidered dynamical range for the anti-Stokes side. The 
data analysis was done using ILL software tools and the 
dynamical structure factor S(Q,E), diffractograms and 
the generalized density of states (GDOS) were evaluated 
using standard ILL procedures without applying multi- 
phonon corrections. 

We complement the experiment with density func- 
tional based ab-initio lattice dynamical calculations to 
(i) reproduce the observed generalized phonon density of 
states (GDOS) and (ii) to extract the partial (atomic) vi- 
brational components contributing to the total spectrum. 
The starting geometry for the calculations was the exper- 
imentally refined structure in the rhombohedral phase 
R3c space group. A G-type magnetic ordering was as- 
sumed as a spin configuration to account for magnetism 
of Fe cations. Relaxed geometries, total energies and 
phonon frequencies were obtained within the K first 
principles framework (which allows a correct treatment 
of the ground state), using similar computational pro- 
cedure as described previously 2 ?. In order to determine 
accurately all the force constants, the supercell approach 
was used for lattice dynamics calculations. A supercell 
(8xa, 8xb, 8xc) was constructed from the relaxed geom- 
etry containing 16 formula units (80 atoms). Total ener- 
gies and HcllmannFeynman forces were calculated for 18 
structures resulting from individual displacements of the 
symmetry inequivalent atoms in the supercell, along with 
the inequivalent cartesian directions ±x, ±y, and ±z. 

The temperature evolution 2 ^ of the dynamic structure 
factor S(Q, E) (Figured]) reveals considerable changes up 
to 1200 K. As a first step it is useful to concentrate on 
the elastic response in order to verify the phase proper- 
ties. The diffraction patterns obtained from the time-of- 
flight data are reported on Figure [3(a) • As indicated by 
the complete disappearance of the magnetic Bragg peak 
at 1.4 A -1 the system loses magnetic ordering close to 
650 K. As the temperature increases further the intensi- 
ties of the two Bragg features located at 1.6 and 2.3 A -1 
decrease steadily reflecting the increased thermal motion 
of the ions. At T c (1050 - 1100 K) the diffraction pattern 
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changes dramatically indicating that a structural phase 
transition took place. As the system is heated further 
there is a second abrupt change of the diffraction pattern 
at about 1200 K, i.e. slightly below the decomposition 
temperature Td = 1240 K. These findings are consistent 
with recent powder neutron diffraction studies 2 — report- 
ing the probable existence of three phases a (T Q < Tc), 
P (T C <T /3 <T 7 ) and 7 (T><T 7 <Tb). 

In view of the strong debate concerning the existence 
of these phases in the bulk material BFO, Figure [Hh) 
shows the phonon spectra, in terms of generalized density 
of states (GDOS) over the explored temperature ranged. 
In order to analyze the temperature-dependent dynamic 
response contained in the measured GDOS we use the 
ab-initio simulated spectrum as a reference point. Fig- 
ure (He) compares calculated and observed phonon spec- 
tra in the ferroelectric a phase. In order to compare 
with experimental data, the calculated GDOS was deter- 
mined as the sum of the partial vibrational densities of 
states gi weighted by the atomic scattering cross sections 
(as seen by neutrons) and masses: GDOS=J2i( (J i/Mi)gi, 
where (oi/Mj = 0.264 (O), 0.207 (Fe), 0.044 (Bi); i={0, 
Fe, Bi}). All the observed features are well reproduced 
by the calculations. The simulated phonon frequencies 
and intensity ratios between the different peaks are in 
good agreement with the measurements. Based on this 
we can consider that the experimental data validate the 
computational approach. We notice that our calculated 
GDOS matches also other available theoretical work 2 ^. A 
more detailed look at the calculated phonons of BFO is 
shown in Figure^d). Therein the partial phonon density 
of states (PDOS) are represented. There is a contribu- 
tion from all the atoms in the range 0-14 meV with 
an increasing amplitude as the atomic mass increases, 
i.e, Bi has the highest intensity in this region. The 
low-frequency region is separated from the rest of the 
spectrum by a gap (14-17 meV). Beyond 14 meV the Bi 
cations do not exhibit any dynamics. The Fe and O den- 
sity of states are both strongly present in the mid-part of 
the spectrum spreading from 14 <E< 45 meV. The oxy- 
gens stretch modes dominate at the highest frequency 
end of the spectrum (45 <E< 70 meV). It should be 
noted that available high-temperature optical spectr 
shows a similar trend comparing to the presently reported 
phonon spectra in terms of temperature-dependent dy- 
namical evolution of BiFeOa. 

The phonon dynamics as a function of temperature 
consistently reflects the suite of phase transitions de- 
duced from the diffraction data. Apart from some in- 
dication for a broadening of the spectra in the oxygen 
dominated region beyond 45 meV the magnetic transi- 
tion leaves a rather weak trace in the phonon spectrum. 
It is mainly characterized by the appearance of strong 
paramagnetic scattering at the location of the magnetic 
Bragg peaks when the transition is passed (Figure [T]) . 
The broadening of the phonon bands and concomitant 
filling of the gap continuos when heating within the po- 
lar ferroelectric a-phase (770 - 1100 K). This indicates 



a very strong anharmonicity of the phonon dynamics in 
this phase. The phonon modes are most affected by the 
transition at T c (1050 - 1100 K) where BFO changes 
from the polar ferroelectric a-phase to a paraelectric in- 
termediate /3-phase. This holds in particular for the low- 
frequency part of the spectrum and makes the statement 
that these modes could provide coupling of magnetism 
and ferroelectricity all the more plausible. In the /3 phase 
(1100 - 1170 K), the intermediate frequency part is sub- 
ject to a considerable downshift resulting in a complete 
filling of the gap that was separating the intermediate 
and low-energy regions. The system undergoes a second 
phase transition around 1200 K to reach a third struc- 
tural phase (7) as clearly reflected in the phonon spec- 
tra. The weight of the low-energy phonon modes dimin- 
ishes when compared to the a and f3 phases in favor of 
two new strong features at 20 and 35 meV. These phase 
transitions unravel changes across various ferroic order 
parameters that are separated over a temperature range. 
As ferroelectricity and magnetism in BFO are driven by 
the Bi and Fe sites, respectively, the partial PDOS re- 
veal clearly the frequency ranges where they contribute. 
Interestingly, the ferroelectic-to-paraelectric phase tran- 
sition involves large changes in the dynamics of the Bis- 
muth cations in terms of vanishment of the related off- 
centered displacements. The abrupt change observed in 
both phase transitions (a — > (3 and j3 — > 7) confirm that 
they are of a first order nature. These results provide 
a reliable support to other works using different tech- 
nique a 11 ! 24 ' 25 . 

As already mentioned, the oxygen contribution to 
the observed phonon spectra changes considerably upon 
heating. It is fully smeared out and concomitantly lost 
intensity in the 7-phase. This is indicative of a weaken- 
ing of Fe-0 bonding leading finally to bond breaking with 
the consequence of oxygen release and vacancy creation. 
This result is quite interesting since it provides a strong 
argument for the occurrence of the high-spin to low-spin 
transition mentioned above, and which is a consequence 
of the reduction of Fe 3+ to Fe 2+ to maintain the charge 
balance due to the presence of oxygen vacancies 3 ^ 24 ' 25 . 
Similar dynamic effects have been seen in superionic oxy- 
gen conductors with an open oxygen network structure 30 . 
The newly visible vibrational features in the 7 phase at 
20 and 35 meV could be due to an increase in cell volume 
after Fe reduction (Fe 3+ — > Fe 2+ ) which results in a red 
shift of closely similar modes observed in the a phase at 
25 and 40 meV, respectively. 

In summary we have carried out inelastic neutron scat- 
tering measurements of the phonon spectra of BiFeOa 
which is a unique model system and multiferroic material 
with multiple order parameters. We have found a strong 
temperature evolution that indicates both strong anhar- 
monicity and subtle dependence of phonons on structural 
details. The phonon modes are most affected by the 
transition at the Curie temperature Tc- Three phases 
a, /3 and 7, in agreement with and in support of other 
works 24 * 2 ^, are highlighted by their clearly different vi- 
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brational signatures. Ab-initio calculations provide a 
phonon density of states in good agreement with the 
observations. The atomistic vibrational understanding 
extracted from these calculations allows gaining reliable 
insights into the dynamic behavior of BFO upon tem- 
perature variation. On the electronic side the occurrence 
of a high-temperature metal-insulator and Fe 3+ — > Fe 2+ 
spin transitions, previously inferred via diffraction exper- 
iments^^, looks very plausible and the present results 



are consistent with these findings. The high temperature 
neutron diffractio n 24 ' 25 and optical data 1 ^ have been as- 
sociated with spin changes and also decomposition and 
our results reveal large changes in the phonon density of 
states, consistent with these findings. We estimate that 
these results are of considerable relevance for the field 
of multiferroism in general and for multiferroism in this 
specific material, in particular. They provide valuable ar- 
guments to better interpret and support results obtained 
with other techniques^-.^ 1 ^^^. 
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FIG. 1. (Color online) Temperature-dependent evolution of the experimental Bose-factor corrected dynamical structure factor 
S(Q, E) of BFO. For clarity, a logarithmic representation is used for the intensities; dark red and dark blue refers to strong 
and weak amplitudes, respectively. 
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FIG. 2. (Color online) Temperature dependence of diffraction patterns (a) and phonon generalized density of states (GDOS) 
spectra (b) of BFO extracted from time-of-flight measurements on IN6. Simulated total GDOS and the related partial atomic 
contributions from ab-initio calculations are shown in (c) and (d), respectively. On the same plot (c) the calculated GDOS 
is compared with a measured spectrum. For the sake of clarity diffractograms and all spectra in (a), and (b), (c) and (d), 
respectively, are vertically shifted with respect to each other. For a better visibility, stability regions and characteristic 
temperatures Tjv and Tc are highlighted with different colors on (a) and (b). 



